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ABSTRACT 

We study the effects of Active Galactic Nuclei (AGN) feedback on the formation and 
evolution of galaxies by using a combination of a cosmological A'^body simulation of 
the concordance A cold dark matter (ACDM) paradigm and a semi-analytic model 
of gal axy formation. This model is an improved version of the one described bv ICoral 
(|2006D . which now considers the growth of black holes (BHs) as driven by (i) gas 



accretion during merger-driven starbursts and mergers with other BHs, (ii) accretion 
during starbursts triggered by disc instabilities, and (iii) accretion of gas cooled from 
quasi-hydrostatic hot gas haloes. It is assumed that feedback from AGN operates 
in the later case. The model has been calibrated in order to reproduce observational 
correlations between BH mass and mass, velocity dispersion, and absolute magnitudes 
of the galaxy bulge. AGN feedback has a strong impact on reducing or even suppressing 
gas cooling, an effect that becomes important at lower redshifts. This phenomenon 
helps to reproduce the observed galaxy luminosity function (LF) in the optical and near 
IR bands at z = 0, and the cosmic star formation rate and stellar mass functions over a 
wide redshift range (0 < z < 5). It also allows to have a population of massive galaxies 
already in place at z > 1, which are mostly early- type and have older and redder 
stellar populations than lower mass galaxies, reproducing the observed bimodality 
in the galaxy colour distribution, and the morphological fractions. The evolution of 
the optical QSO LF is also reproduced, provided that the presence of a significant 
fraction of obscured QSOs is assumed. We explore the effects of AGN feedback during 
starbursts finding that, in order to obtain a good agreement with observations, these 
need to be strong enough to expell the reheated gas away from the galaxy halo. 
We also test new, recent prescriptions for dynamical friction time-scales, and find 
that they produce an earlier formation of elliptical galaxies, and a larger amount of 
disc instabilites, which compensate the change in the merger frequency such that the 
properties of z = galaxies remain almost unaffected. 
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general 



1 INTRODUCTION 

Over the last decade, there has been a phenomenal increase 
in the range and quantity of data on quasi-stellar objects 
(QSOs) or 'quasars' and galaxies with active galactic nu- 
clei (AGN), which are the seats of an intense production of 
energy. They were singled out from the general galaxy pop- 
ulation by their peculiar colours, morphology and/or strong 
emission in the radio. X-ray or infrared bands. The conven- 
tional view regards QSOs and AGN as different manifesta- 
tions of the same phenomenon. Observations of the brightest 
Seyfert galaxies, one of the types in which AGN are classi- 



fied, show events of sudden bursts of star formation activ- 
ity (starbursts) that occurred 1-2 Gyrs ago; these events 
have been associated, by some authors, with galaxy mergers 
(|Kaufrmannll2003l; ISanders fc Mirabeilll996l ). Recent results 
by Feain et al. ( 20071 ) show direct evidence of the QSO-star 
formation connection where 70% of the quasar radio emis- 
sion comes from star formation activity of the host galaxy. 

In a pioneering work bv lHovle fc Fowled (ll963D . it was 
proposed that the energy from AGN was of gravitational ori- 
gin derived from the collapse of very massive objects under 
their own strong gravitational fields. These early ideas found 
a modified expression in the paradigm of black hole accretion 
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discs a few years later l|Lvnden-Belllll969l ). In more recent 
times, an important connection between AGN and galaxy 
evolution has been deduced from the correlation between 
black holes, thought to be the engines of AG N, and proper- 
ties o f the host galaxy such as the bulge mass (|Haring fc Rixl 
I2OO4I ), the luminosity iii the K -^ B- and V^-bands, the stellar 
mass arconi fc Huntll2003l). and the velocity dis persion of 
the stars in the bulge (jFerrarese fc Merritll2OO0l ). The ex- 
istence of a galactic bulge has been thought to be related 
to merger events, which led to speculate that mergers may 
drive the formation and evolution of black holes. Further- 
more, iBovle fc TerlevichI { 1998h show a correlation between 
the evolution of the global stellar formation rate and the 
luminosity density of optically selected quasars. This could 
also indicate that black holes are related to merger events 
and, therefore, to s tarburst phe nomena as many authors 
have speculated (e .g. Malbon et al . 2007; Croton ct al. 2006; 
iBower et al]l2006l : ISiiacki et alH?007l 1 . 

The alleged link between AGN and star formation ac- 
tivity brings questions related to the effect of this nuclear 
activity on the evolution of their host galaxies or galaxy clus- 
ters. Many recent observational results have began to answer 
some of these questions complementing our picture of galaxy 
forma tion and evolution. For example, ISchawinski et all 
iiooj) find a population of early type galaxies lying outside 
the red sequence in the colour-magnitude diagram; this red 
population is comprised by star forming and Seyfert galax- 
ies. From these data, they identify an evolutionary sequence 
from star forming (blue population) to quiescent galaxies 
(red population) . This transition is explained by the effect of 
an AGN phase occurring roughly 0.5 Gyr after the star burst, 
which suppresses star formation. iReuland et"al] (|2007l ) have 
been able to confirm this result at high redshifts. 

Direct observational evidence of AGN activity is pro- 
vided by radio galaxies and quasars containing jets that 
transport energy from the AGN to the surrounding medium. 
The effects of jets issuing from AGN have been observed 
in many different wavelengths. Radio sources show large- 
scale outflows of si gnificant gas masses which are gener- 
ated by radio jets (INesvadba et al.l I2OO6I : iNesvadba et al.l 
l2007l : lTemi et al.l [20071 ). Theoretical studies of jets indicate 
that the time-scales for transient activity (e.g. restoration 
of equilibrium, buoyant transport in the hot gas) are consis- 
tent with a release of feedback energy from a central black 
hole (|Temi et al.ll2007l ). The dichotomy in the properties of 
low- and high-excitation radio galaxies has been interpreted 
as an additional contribution from cold gas acting in high- 
excitation source s that forms the co ld disc and torus around 
a central engine (|Evans et al.ll2007l '). Observations in X-rays 
indicate the presence of shocks, bubbles and sound waves 
thought to be driven by nuclear activity and to be respon- 
sible for heatin g the cluster core, as in the Perseus cluster 
l|Sandersll2007l . and references therein). 

All these results illustrate the key role of active 
galactic nuclei on galaxy evolution. As AGN have be- 
come a crucial ingredient of the galaxy evolution pro- 
cess, many theoretical works have proposed different pos- 
sible scenarios for BH growth and the associated AGN 
feedback. These works include pure analytic approxima- 
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zuua ; >_.roron er ai.i i2006 ; Bower et al.l |2006 : , Malbon et al 
2OO7I : iMaruUi et all [20081). and fu l ly self-consistent n umer 



ical models (jSpringel et al.ll2005al : [siiacki et al.ll2007l ). The 

application of a semi-analytic model of galaxy formation on 
a cosmological numerical simulation has demonstrated to 
be a very appropriate tool for the study of phenomena that 
drive the formation and evolution of galaxies within the hi- 
erarchical clustering scenario. This is the method used to 
investigate the effects of AGN feedback in the present work. 

Our semi-a nalytic model is based on the one described 
bv ICoral i 20061 ) ■ to which we refer to as SAGl (acronym 
for 'Semi- Analytic Galaxies' version 1). This model follows 
the formation and evolution of galaxies including gas cool- 
ing, star formation, SNe explosions and mergers, as well 
as a detailed implementation of metal enrichment of stars 
and interstellar and intergalactic medium. SAGl, as well 
as many other semi-analytic models, reproduce to a rea- 
sonable degree observational results such as gas fractions, 
luminosity functions (LF), the TuUy-Fisher relation, and 
galaxy colours. However, some inconsistencies still persist, 
in particular, the existence of a blue massive galaxy popu- 
lation and an excess of bright galaxies in th e LF compared 
with observations in the fej- and K-hand ( e.g. lSpringel et al.l 
I2OOII : ICole et al.ll2000l : iBaugh et all 12004 ). These problems 
are in part due to the fact that the mass function of host 
dark-matter haloes does not show a sharp cut-off at masses 
corresponding to L* galaxies in the model. Therefore, mod- 
els that popula te these halo es tend to produce too many 
bright galaxies (|Baughll2006l ). The implementation of an ar- 
tificial supression of gas cooling made in SAGl, that affects 
haloes characterized by a virial velocity higher than a cer- 
tain threshold, was a simple solution to overcome, at least 
partially, the inconsistencies mentioned above. These prob- 
lems call for the inclusion of additional processes in order 
to reconcile, to a further degree, model results with obser- 
vations. The physical process thought to be responsible for 
the control of the amount of gas that can cool is AGN feed- 
back, which is triggered by gas accretion events onto central 
supermassive black holes. Recent works have made progress 
in this sense, demonstrating that AGN feedback is crucial to 
reduce gas cooling in large haloes, thus preventing them from 
forming stars at late times and being able to explain the ex- 
ponential cut-off in the the bright-end o f the galaxy LF (e.g . 



Cattaneo et"alll2006l: iMenci et al.|[200^ : ICroton et al.ll2006l : 
Bower et al.ll2006l : iMalbon et al.ll2007l ). Following the pre- 



tions (e.g. lEfstathiou fc Reeslll988l: iPercival fc Milleijll999l : 
[Ch urazov et al. I I2OO2I : iHopkins et all |2007| ), semi-analytic 
models of galaxy formation 



scriptions implemented in these works, we replace the simple 
procedure of stopping gas cooling in haloes with high virial 
velocity, as it is implemented in SAGl, by a more physical 
modelling of AGN feedback in galaxies, giving raise to the 
version SAG2. 

The currently implemented models of AGN feedback 
only prevents significant gas cooling in large haloes. The so 
called radio mode feedback (see for instance, ICroton et al.l 
l200d) is assumed to produce low energ y 'radio ' activ ity. Re- 
cent observational results bv G astadello et al.l l|2007h on the 
galaxy cluster AWM 4 cannot be explained if AGN activ- 
ity only produces this mode of feedback, since it presents 
an active radio source which is not accompanied by a X-ray 
temperature gradient; this might indicate the presence of 
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AGN feedback during high accretion rates originated dur- 
ing mergers and disc instabilities, which also feed the BH 
growth. We explore this possibility within the A cold dark 
matter (ACDM) framework using SAG2. 

Our model SAG2 is also used to investigate the impact 
of new theoretical merger t ime-scales p r esente d in two par- 
allel studies carried out bv I Jiang et alj ()2007h and Boylan- 
Kolchin, Ma & Quataert (2007). These wo rks show that 
the eq uation for merger time-scales given bv lLacev fc Colj 
(| 19931 ) systematically underestimates the time-scales of mi- 
nor mergers and overestimates those of major mergers. 

This work is organised as follows. Section [2] presents 
a description of the cosmological simulation used and the 
improved semi- analytic model SAG2, summarising first the 
original model SAGl and then giving a detailed description 
of the AGN feedback model introduced. Section [3] describes 
the behaviour of BH and QSO properties in the model, 
and presents the analysis of the effects of AGN feedback 
on galaxy properties. Further improvements made to SAG2, 
including the effects of starburst AGN feedback, and new 
prescriptions for merger time-scales, are presented in Sec- 
tion |4l Finally, Section [5] gives the main conclusions of this 
work. 



2.1 ACDM Cosmological simulation 

We use a cosmological simulation of the concordance ACDM 
cosmology in a periodic box of 60 Mpc. This box is 
large enough to allow us to make a suitable analysis of 
the luminosity function of galaxies and quasars, and of the 
galaxy mass function. The simulation contains 16, 777, 216 
dark matter particles with a mass resolution of 1.001 x 
10^ h~'^MQ. More than 54,000 dark matter haloes have 
been identified, with the largest one having a mass of 
5.36 X 10" h-^ Mq. 

The simulation p arameters are consistent with the re- 
suhs of WMAP data (|Spergelll2003h . that is, fim = ^DM + 
f^baryons ~ 0.28 (with a baryou fraction of 0.16), flA ~ 0.72 
and (T8 = 0.9. The Hubble constant is Ho = 100/iMpc"\ 
with h — 0.72. The gravitational softening length, e, is 
3.0 kpc. The simulation starts from a redshift z = 48 
and is run using the public version of the GADGET-2 code 
(|Springelll2005b[ ). 

Dark haloes are first identified as virialized particle 
groups by a frien ds-of-friends (FQF) algorithm. A SUB- 
FIND algorithm jSpringel et al.1 bOOll ) is then apphed to 
these groups in order to find self-bound dark matter sub- 
structures. Merger trees are then constructed from these 
dark matter haloes and their embedded substructures. 



2 A HYBRID MODEL OF GALAXY 
FORMATION AND EVOLUTION 

We study the formation and evolution of galaxies by apply- 
ing a numerical technique which combines a cosmological 
iV-Body simulation of the concordance ACDM universe and 
a semi-analytic model of galaxy formation. In this hybrid 
model, the outputs of the cosmological simulation are used 
to construct merger histories of dark matter haloes which 
are used by the semi-analytic code to generate the galaxy 
population. The main advantage of semi-analytic codes is 
that they allow to reach a larger dynamic range than fully 
self-consistent simulations, at a far smaller computational 
cost. 

The s emi-analytic model used here is based on that de- 
scribed bv lCoral (|2006l ) (SAGl). The physical processes con- 
sidered in this model are cooling of hot gas as a result of 
radiative losses, star formation, feedback fro m supernova ex- 
plosio ns and galaxy mergers, modelled as in ISpringel et all 
l|200ll ). The model also tracks the circulation of metals be- 
tween the different baryonic components, that is, cold gas, 
hot diffuse gas, and s tars, following the model described by 
iDe Lucia etal] l|2004l ). with the additional advantage of in- 
cluding the mass evolution of different chemical elements. 

In the present work, SAGl is improved by including 
feedback from AGN as a replacement of the artificial sup- 
pression of gas cooling in massive haloes, thus leading to 
the version SAG2. In the following subsections, we present 
the main characteristics of the iV-Body simulation used, and 
briefiy describe the physical processes already included in 
SAGl. In the original version of the code, only major merg- 
ers trigger starbursts; we now include two other mechanisms 
responsible of this kind of events, minor mergers and disc in- 
stabilities. Then, we present the model of black hole growth 
and the associated AGN feedback implemented in the new 
version of the code. 



2.2 Semi-analytic model SAGl 

We briefiy comment the processes of gas cooling, star for- 
mation, supernova feedback and metal pro duction inclu ded 
in the model SAGl; we refer the reader to ICoral (|2006l ) for 
a more detailed description. 

We assume that hot gas always fills the dark matter 
haloes following an isothermal distribution. Initially, the hot 
gas mass is given by the baryon fraction of the virial mass 
of the dark matter halo. It is subsequently modified by the 
amount of gas that has cooled and the stellar mass formed 
in the galaxies contained within the halo. The mass of hot 
gas that cools in each dark matter halo in the underlying 
cosmological simulation is given by the cooling rate 



dMcool . 2 d^cool 

= 47rpgrcooi - 



dt 



dt 



(1) 



where pg is the density profile of an isothermal sphere, and 
J'cooi is the cooling radius. The local cooling time is defined 
as the ratio between the specific thermal energy content of 
the gas and the cooling rat e per unit volume, A(r, Z), which 
depends of the metallicity l|Sutherland fc Dopita|[l993 ) , and 
the temperature T = 35.9(yvir/kms~^) of the halo, being 
Vvir its virial velocity. This process only operates when the 
galaxy is the central galaxy of a given halo. The cooled gas 
contributes to the formation of the galaxy disc. The cold 
gas mass of each galaxy is involved in the star formation 
process. The star formation rate is given by. 



(2) 



dt ^ tT' 

dyn 

where a is a parameter which regulates the efficiency of 
star formation, McoidGas is the cold gas mass, and = 
0.1i?vir/Vvir is the dyuamical time of the galaxy, being Rvir 
the virial radius of the halo. Here a = ao(K,ir/220kms~^)", 
where qq and n are free parameters set to 0.1 and 2.2, 
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respectively. Each star formation event generates a stellar 
mass AAf*, which leads to a number rjcc of core collapse 
supernovae (SNe CC). This quantity depends on the ini- 
tial mass function (IMF) adopted, which in this case is a 
Salpeter IMF normalized between 0.1 and 100 M©; there- 
fore, 7?cc = 6.3 X 10"^ Mg^ SNe CC include types Ib/c 
and II. The energy Esncc released by each SN CC (1.2 x 
10^^ ergs~^) is assumed to reheat the cold gas of a galaxy 
inducing galactic outflows that transfer cold gas to the hot 
phase by the amount 



now implemented in the semi-analytic code as we describe 
in Section \2A\ 



AAf _ 4 ?7cc-EsNCC . 

AAJrchcatcd — TTj AM* 



(3) 



where e is a dimensionless parameter which regulates the 
efficiency of the feedback by SNe CC, and takes a value 
of e = 0.1. The reheated gas is k ept within its host d ark 
matter halo ('retention model', see lDe Lucia et aLll2004l for 
other possibilities). 

We set an initial primordial abundance of the baryonic 
components consisting of 76% hydrogen and 24% helium. 
They become chemically enriched as star formation and 
metal production take place. Stars contaminate the cold and 
hot gas as a result of mass losses during their evolution. A 
fraction /ojec of the mass ejected by galaxies is transferred 
directly to the hot phase instead of being first incorporated 
to the cold gas. We adopt here /cjoc ~ 0.5, a value that allows 
to obtain an evolution of the ma ss-metallicity re lation con- 
sistent with observational data (|Erb et al.ll2006l ). resulting 
in a more gradual contamination of the cold gas associated 
to each galaxy. The transfer of reheated enriched cold gas 
also contributes to the chemical enrichment of the hot phase. 

Since the cooling rate depends on the metallicity, this 
chemical enrichment has a strong influence on the amount 
of hot gas that can cool. This process in turn influences 
the star formation activity which is ultimately responsible 
for the chemical pollution. Thus, new stars formed from the 
polluted cold gas become more chemically enriched. Stars 
not only form in a quiescent way from the reservoir of cold 
gas (Eq. [2} but also from starbursts triggered by mergers 
and disc instabilities, as we will describe in Section [2]3] The 
chemical model implemented in our semi-analytic code con- 
siders that metals are provided by three kinds of sources: 
low and intermediate mass stars, SNe CC, and type la su- 
pernovae (SNe la); the flrst group of sources yields metals 
through m ass losses and stellar winds (details are given in 
ICorall2006h . 

With respect to the spectro-photometric properties 
of galaxies, we consider evolutionary synthesis models 
that depend on the metallicity of the cold gas from 
which stars are formed, estimated for a Salpeter IMF 
l|Bruzual fc CharlotI 12003 ). All magnitudes and colours in- 
clude the effects of dust extinction following the implemen- 
tation bv IPe Lucia ct al. (20041. which is based on previous 
works (e.g. Kauffmann et al.lll999al ). 

It is important to remark that SAGl suppresses gas 
cooling in rapidly rotating dark haloes (Vvir > 350kms~^). 
This restriction is set with the aim of avoiding very mas- 
sive and luminous galaxies, thus achieving a better agree- 
ment with observations of the bright end of the luminosity 
function. The reduction or suppression of th e cooling flow 
can be naturally produc ed by AGN feedback (jCroton et al.l 
l2006l . lBower et aLll2006l ). This important physical process is 



2.3 Starbursts in the model: Galaxy mergers and 
Disc Instabilities 

Our new version of the semi-analytic model produces star- 
bursts in three different ways: major mergers, minor mergers 
and disc instabilities. Starbursts occurring in the last two 
processes are new ingredients added to the original version 
SAGl. All the cold gas contained in the galaxies involved 
in any of these three processes is completely consumed in a 
sudden burst of star formation. 

In a hierarchical scenario of structure formation, merg- 
ers of galaxies are a natural consequence of the mergers of 
dark matter haloes in which they reside, and play an im- 
portant role in determining the mass and morphology of 
galaxies. In the subhalo scheme arising from the identifica- 
tion of dark matter substructures, we distinguish three types 
of galaxies when tracking galaxy formation. The largest sub- 
halo in a EOF group hosts the 'central galaxy' of the group; 
its position is given by the most bound particle in that 
subhalo. Central galaxies of other smaller subhaloes con- 
tained in a EOF group are referred to as 'halo galaxies'. 
The (sub) haloes of these satellite galaxies are still intact af- 
ter falling into larger structures. The third group of galaxies 
comprises 'satellite' galaxies that arise when two subhaloes 
merge and the galaxy of the smaller one becomes a satellite 
of the remnant subhalo. In SAGl, satellite galaxies are as- 
sumed to merge on a dynamical time-scale with the central 
galaxy of the largest subhalo of the FOE group they reside 
in. It is assumed that the satellite galaxy has a circular or- 
bit with a velocity given by the virial velocity of the parent 
halo and decays to the corresponding central g alaxy through 
dynamical friction l|Binnev fc TremainelllQSTf) , given by 



Ttric 



1/W 



(4) 



2 C GM''^'' InA' 

where Mgat is the satellite mass orbiting at a radius rc inside 
an isothermal halo of circular velocity Vc. The function /(e) 
describes t he dependence on th e eccentricity of the satellite 
orbit (see ISpringel et al] l200ll for a detailed description), 
and C = 0.43. The function InA = {1 + m'^™"^'''/^^'''* ) is 
the Coulomb logarithm. 

Mergers are classified according to the ratio between 
the mass of the accreted satellite galaxy and the mass of 



the central galaxy, /merge = M /M" 



A major merger 



occurs when /merge > 0.3; in this case all stars present are re- 
arranged into a spheroid and all the cold gas in the merging 
galaxies is assumed to undergo a starburst. The stars thus 
produced are added to the spheroid (bulge). Major mergers 
constitute a mechanism to produce elliptical galaxies. Sub- 
sequent accretion of cooling gas is tunneled into the central 
galaxy generating a disc. 

On the other hand, the merger is classified as a mi- 
nor merger when one of the two galaxies has a small mass 
compared to the o ther, that is, /merge < 0.3. Following 
iMalbon et al.l l|2007l ). we assume that the presence of a star- 
burst in a minor merger depends on the gas mass fraction of 
the disc of the central galaxy, fg'Zal = MS-J^SL/MdTJL*^"', 
where Mdisc = Msteiiar - MBuige + McoidGas • In this last rela- 
tion, Afsteiiar IS the total stellar mass of the galaxy, Afsuige is 
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the mass of the bulge (formed only by stars) , and McoidGas 
is the cold gas content of the galaxy. If /coWgLs > /gas.burst , 
where /gas.burst is a free parameter of the model, the mi- 
nor merger will produce a starburst as a result of the per- 
turbation introduced by the merging satellite galaxy which 
drives all the cold gas from both galaxies into the spheroid, 
where it is completely transformed into stars. Otherwise, 
if /coMgL < /gas.burst then no burst occurs. We adopt a 
threshold /gas.burst = 0.6. Furthermore, if the satellite is 
much less massive than t he central galaxy (/m erge < ./burst, 
with /burst = 0.05, as in iMalbon et al.ll2007l ). the central 
stellar disc remains unchanged and the accreted stars are 
added to the spheroid; in this case, there is no starburst, 
regardless of the amount of gas in the central disc, which 
already contains the contribution of the gas of the satellite. 

Starbursts are also triggered by disc instabilities. When 
a galaxy disc is sufficiently massive that its self-gravity is 
dominant, it becomes unstable to small perturbations by 
minor satelli tes or dark mat ter substructures (Mo, Mao & 
White 19 98:ICole et al.ll2000l ). We follow the stability crite- 
rion from lCole et al.l "( 20001 ). given by 



(GMdisc/rdisc)i''2 ' 



(5) 



where Vmax is the maximum circular velocity of the disc, 
and Mdisc is the disc mass. We assume that raise is the disc- 
scale radius (following IMo et al.lll998l ). When e < tdisc the 
disc becomes unstable; in the presence of a perturbation 
all the stars and cold gas in the disc are transferred to the 
bulge. In this case, all the gas present in this modified bulge 
is assumed to be consumed in a starburst. After this pro- 
cess, discs are required to stabilize before a new starburst 
by disc instability can be triggered. We adopt edise = 1-1 
taking into account different observational contraints, as we 
describe subseque ntly along t he pa per. This value lies within 
the limits given bv lMo et all (Il998l ') (ranging from 0.8 to 1.2) 
for dark matter haloes with NFW profiles. This mechanism 
also generates elliptical galaxies. 

The upper panel of Fig.[T]shows the evolution with red- 
shift of the number density per unit volume of events that 
produce starbursts, including major mergers, minor mergers 
with starbursts, and disc instabilities. As can be seen, star- 
bursts arising from major and minor mergers peak around 
~ 2 and 3, respectively. Disc instability events increase 
in number as the redshift decreases, with a number den- 
sity that peaks at z ^ 1. At lower redshifts, the number of 
disc instabilities is higher than that of minor mergers. The 
lower panels of Fig. [T]show the number of events that trigger 
starbursts weighted by the ratio between the mass of stars 
produced during the burst and the stellar mass, for galaxies 
of high stellar meiss, M* > 10^° Mq (lower left panel), and 
low stellar mass, M^, < 10^*^ M© (lower right panel). It can 
be seen that, regardless of the stellar mass, disc instabilities 
are the preferred mechanism for the growth of stellar mass 
at all redshifts. Major and minor mergers play a minor role 
in the growth history of the stellar mass. An important point 
to notice is that the growth of stellar mass in high mass sys- 
tems reaches a peak at redshifts z ~ 1.5, whereas low stellar 
mass galaxies continue to grow almost down to 2: = 0. 
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Figure 1. Upper panel: Evolution with redshift of the number 
density of events that trigger starbursts in our model; major merg- 
ers (dashed lines), minor mergers with /gas > /gas.burst (dotted 
lines), and disc instabilities (solid lines). Lower panels: Evolution 
with redshift of the number of events weighted by the ratio be- 
tween the mass of stars produced by the burst and the stellar 
mass present at the moment of the burst for high (lower left) and 
low (lower right) stellar mass galaxies. Lines are as in the upper 
panel. 



2.4 Black Hole growth, AGN outflows and gas 
cooling suppression 

As it was discussed before, there is growing evidence that 
AGN are extremely important for galaxy formation and evo- 
lution. AGN are candidates to suppress or reduce cooling 
flows, thus afi'ecting the star formation process and asso- 
ciated galaxy properties, such as the luminosity function, 
colours, and ste l lar m asses. As it was first pointed out by 
iHovle fc Fo^ l| 19631 ). the energy source of a quasar or 
AGN is gravitational and could arise from a highly collapsed 
object or a massive black hole, an idea broadly agreed to by 
most authors in the field. Thus, in order to produce AGN, we 
first need to follow the growth of central black holes within 
galaxies, for which we propose a physi cal model based on a 
co mbination of prescript ions given by ICroton et al.l (I2OO6I ) 
and lMalbon et all |2003). Considering that QSOs and AGN 
are different manifestations of the same phenomenon, the 
models of BH growth used in these works are motivated 



6 Claudia del P. Lagos, Sofia A. Cora, Nelson D. Padilla 



by the observational evidence linking the global star forma- 
tion rate and the optica l QSO density in the universe (e.g. 
iBovle fc Terlevichlll998l ) in which both are correlated. Given 
that there are two distinct ways to form stars, either quies- 
cently or in bursts, we also make a difference between the 
way in which a BH grows in these two modes. We describe 
these modes separately in the following two subsections. 



2.4-1 BH growth in violent star formation processes 

We assume that a BH grows from perturbations to the 
gaseous disc resulting from galaxy mergers or disc insta- 
bilities. In such cases, cold gas is funncl ed to the cen tre of 
the galaxy producing starbursts l|Croton et al.l [20061 '). The 
growth of BIfs as a result of these processes is referred to as 
'starburst mode growth'. During a galaxy merger, we assume 
that the central BHs merge instantaneously. For simplicity, 
we ignore the effects of gravitational waves in BH mergers 
so that the resulting mass is the direct sum of the two merg- 
ing masses, with the BH mass in the satellite galaxy being 
transferred to the BH mass in the central one. The resulting 
object continues to accrete cold gas. The mass accreted by 
the BH is proportional to the total cold gas mass present; 
the efficiency of accretion is lower for unequal mergers. This 
is summarised by 



AMb 



/bh 



A^ColdGa 



1 + (200kms-VKir 



(6) 



where McoidGas is the cold gas mass of the central galaxy 
after adding the cold gas of the satellite. Note that wc use 
a velocity of 200kms~^ in the above formula instead of a 
value of 280kms~^ as considered by ICroton et al.l l|2006l ). 
As in that work, the value of the parameter /bh is chosen 
so as to reproduce the observed local relation between the 
BH mass and the host bulge mass (Mbh — Mbuigo relation), 
as we discuss in Subsection 13.21 we adopt /bh = 0.015. 

The BH growth through disc instabilities is also repre- 
sented by Eq.ia where the ratio m=''7M'^''°"'^' is replaced by 
unity, since this process is a dynamical self-interaction that 
depends on the properties of a single galaxy. The behaviour 
of the BH growth during starbursts is shown in Fig. [5J and 
will be discussed in Subsection 13.11 

It is important to note that there are no BH seeds in 
our model and, therefore, the birth of central super massive 
BHs is triggered when galaxies undergo their first starbursts 
at some point in their evolution. 



2.4-2 BH growth in gas cooling processes 

The second way in which BHs are assumed to grow is from 
cold gas accretion during gas cooling. This occurs once a 
static hot gas halo has formed around the host galaxy. Note 
that satellite galaxies do not experiment this growth mode 
since the reservoir of hot gas is associated only to the central 
galaxy. 

Following ICroton et al.l (|2006l ) , we assume this gas ac- 
cretion onto the BH to be continuous and quiescent and to 
be described by a simple phenomenological model. 



Mbh = KAGN 



Mbh /hot 
IO8M0 ^ "oT ^ 



Vvir 



200 kms-i 



(7) 



In this equation, Mbh represents the black hole accretion 
rate, AIbh is the black hole mass, /hot is the fraction of the 
total halo mass in the form of hot gas, /hot ~ mnotGas/Mvir, 
where A/vir is t he virial mass of the host halo, and kagn is a 
free parameter. ICroton et al.l l|2006l) note that /hot becomes 
approximately constant when Vvir ^ 150kms~^; the depen- 
dence of A^BH on this quantity is not important. 

The mechanical heating generated by the BH accretion 
(black hol e luminosity) can be expressed, following the ar- 
gument of lSohanI (|1982| '). as 



Lbh = V Mbhc 



(8) 



where c is the speed of light, and rj is the standard efficiency 
of energy production that occurs in the vicinity of the event 
horizon. This process may inject sufficient energy into the 
surrounding medium to regulate gas cooling since it reduces, 
or even stops, the cooling flow. Consequently, the cooling 
rate has to be modifled accordingly. 



M^. 



Afc, 



vir / 



(9) 



This feed back process is referred to as 'radio mode feedback', 
following ICroton et al.l (|2006l ) ; the way in which BHs grow 
during this process is thus called 'radio mode growth'. 

The BH luminosity depends on the combination of the 
free parameters kagn and r\ (Eq.[7]and[Sl respectively). The 
local AfBH — Afbuigc relation does not help in tuning these 
values since it is not sensitive to them. However, the QSO 
and galaxy LFs are good observational constraints. We find 
that by adopting kagn = 2.5 x 10"* M© yr~^ and r\ — 0.1, 
the observed shape of the QSO LF and the bright-end break 
of the galaxy LF are fairly well reproduced as it is described 
in Subsections 13.31 and 13.4.11 respectively. 

Notice that we are not taking into account the fact that 
BHs may be radiatively inefficient when accreting gas during 
the cooling flow process. This possibility is supported by the 
presence of BHs at the centre of clusters with cool cores 
that are usually weak radio sources (e.g. M87 in the Virgo 
Cluster, iLe fc Becked l2007l . IWang et al.ll2008l ). While this 
particular phenomenon might be accounted for by a non- 
constant efficienc y r;, some models use l ower BH accretion 
rates instead fe.g. IJoUev fc Kuncidl2007l . IStern fc PoutanenI 

boosh. 



2.4-3 Remarks on BH luminosity, Lbh 

The definition of the black hole luminosity, Lbh, given in 
the previous section (Eq.[S}, involves the mass accretion rate 
onto the black hole. This quantity is given by Eq. [7] for the 
radio mode BH growth. In the case of the starburst mode 
growth, which occurs during mergers and disc instabilities, 
the mass accretion rate is estimated from the accreted BH 
mass in the starburst mode, as it is given by Eq. |6] 



Mb 



AMb 



AT 



(10) 



where AT is the time interval determined by the 50 subdivi- 
sions between consecutive simulation snapshots, introduced 
to increase the accuracy of the integration of differential 
equations; this is the same time interval applied to estimate 
AAf* from Eq. [3 
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In this way, the luminosity of the black hole will grow 
every time the galaxy undergoes star formation in mergers 
and disc instabilities, and during gas cooling processes with 
soft, quiescent formation. 



2.5 The Eddington Limit in the 'BH' Luminosity 

For a black hole of mass Mbh, the maximum allowed lumi- 
nosity is the Eddington luminosity, defined as 

LEdd = 47rGMBH^, (11) 

where G is the gravitational constant, rup is the proton mass 
and (tt is the Thompson cross-section for the electron; the 
latter assumes a gas composed mostly by hydrogen. In this 
expression, we consider the black hole mass at the end of the 
accretion process. When the calculated black hole luminosity 
is greater than the corresponding L^dd, we set Lbh = L^dd- 
Nonetheless, we do not limit the BH accreted mass, a choice 
justified by the lack of strong observational constraints on 
the efficiency rj, that regulates the amount of accreted mass 
that is transformed into black hole luminosity. In our model, 
the maximum instantaneous fraction of BHs with luminosi- 
ties over the Eddington limit is only ~ 6% at very high 
redshift, so this assumption implies a minor effect. 



3 BEHAVIOUR OF THE NEW MODEL SAG2 

The principal aim of introducing AGN feedback in models 
of galaxy formation is to account for the observed influence 
of AGN on the evolution of galaxies. I n particular, the inclu- 



200e 



200e 



sion of AGN in semi-analytic models (ICattaneo et al.l 
Menci et al.' '2006'; Croton ct al. '2006*: ' Bower et al.l 
Malbon ct al. 2007; Marulli ct al. 2008) has also proven to 
be a good way to avoid problems present in previous ver- 
sions of semi-analytic models such as an excess of bright, 
extremely blue galaxies, when compared to observations (as 
occurs in SAGl). 

In the original model SAGl, a number of free param- 
eters was set to regulate the gas cooling, star formation, 
supernovae feedback and galaxy mergers, and to determine 
the circulation of metals among the different baryonic com- 
ponen ts. These parameters were carefully tuned in ICoral 
l|2006l ) to satisfy numerous observational constraints such 
as Milky Way properties, the luminosity function, and the 
TuUy-Fisher (TF), colour-magnitude and mass-metallicity 
relations. An important success of the original code resides 
in that the implemented chemical model was able to explain 
the spatial distribution of iron in the intracluster medium at 
z = 0, and the te mporal evolution of its mean iron content 
ICora et al.ll2008l ). 

Once the new free parameters involved in the black hole 
growth and AGN activity are added, a new calibration of the 
model is performed to reproduce the above mentioned ob- 
servational constraints plus new observational relationships 
between the properties of black holes and their host galax- 
ies, that is, BH mass-bulge mass, BH mass-bulge luminosity 
relations, and the QSO luminosity function. 

In the remainder of this section, we analyse the model 
predictions on BH relations and how the heating of AGN 
affects galaxy properties. 



3.1 Black hole growth 

The three different processes driving the BH growth are gas 
cooling, galaxy mergers and disc instabilities. Their influ- 
ence on the volume averaged individual BH growth is shown 
in Fig. [2j represented by different line styles. The left panel 
involves all BH masses, showing that disc instabilities are 
the major contributors to the BH growth at all redshifts, 
closely followed by mergers, both with a smooth variation 
over a wide range of redshifts. Gas cooling processes become 
more important at low redshifts with the growth rate steeply 

increasing by two orders of magnitude for z < 2^ 

These results differ from those reported bv lBower et al.l 

(|2006t ). who find that mergers have a small contribution to 
the BH growth, and that disc instabilit ies are the m o st im - 
portant contributors at high redshifts. iBower et al.l (|2006l ) 
also find that gas cooling becomes dominant at low red- 
shifts, a£d_tlmtJt_opeiute since earlier epochs. On the other 
hand, ICroton et al.l (|2006l ) find that the radio mode growth 
is less important than their 'QSO mode' growth (as sociated 
with s tarbursts) in contrast to the results from iBower et al.l 
(|2006l ). and that it stays constant since z 2, a different be- 
haviour to that seen in our model. These differences demon- 
strate that the model implementation and parametrization 
have a notorious influence on the detailed history of BH 
growth. In spite of this, different semi-analytic models are 
still able to reproduce several observational constraints in- 
volving BH masses, such as the 'BH-bulge' relations de- 
scribed in the next SubsectiorQ. 

Fig. [2] also shows the dependence of BH mass growth 
rate on the three above-mentioned mechanisms for low 
(Mbh < IO^Mq, middle panel) and high (Mbh > 10* M©, 
right panel) BH masses. For low mass black holes, we flnd 
that mergers dominate their growth down to z ~ 2, where 
disc instabilities start to dominate for lower redshifts. The 
effects of gas cooling processes are negligible. In the case 
of high BH masses, the behaviour of the three mechanisms 
remains similar to what is shown by the full BH sample. 

The average increase of BH mass with redshift for dif- 
ferent ranges of BH masses is shown in Fig. [S] High mass 
BHs acquire both, half and 90 per cent of their flnal masses, 
at earlier times than smaller BHs. In other words, mas- 
sive black holes form first. This is in agreement with re- 
cent observational works (e.g. lMarconi et al]|2004l ) that find 
that, by 1.5 < z < 2.5, massive black holes have al- 
ready attained at least one half of their final mass; smaller 
BHs only reach half of their final mass at z < 1. This is 
consisten t with findings on the evolution of the X-ray LF 
of QSOs llCowie. Songaila fc Cohen|[l99^ : lUeda et al.|[2003l : 
iBarger et al.l I2005D . where the number of bright sources 
is higher at higher redshifts. Moreover, observations have 
shown that massiv e galaxies are alrea dy in place in the dis- 
tant universe (e.g. iDrorv et al]|2004 ). With respect to star 
forming galaxies, it has been found that their luminosities 
in the near IR bands decline with decr easing redshift since 
z ~ 1 (|Cowie. Songaila fc Cohenlll996l ). All these observa- 



^ The different BH growth histories presented by these models 
will show important effects in the global galaxy population at 
different redshifts. As we show in the remainder of this paper, 
we use several galaxy statistics to help reduce the degeneracies 
between semi-analytic models. 
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Figure 2. Volume averaged individual BH growth rate as a function of redshift for the three mechanisms that drive the BH growth: 
gas cooling processes (dashed line), galaxy mergers (dotted line) and disc instabilities (solid line), for all BHs masses (left panel), small 
black holes (< lO*Af0, middle panel) and massive black holes (> 10*A/q, right panel). 




Figure 3. Average growth history of BHs for different mass 
ranges starting from 2 = 5. The tracks are constructed for a 
fixed number of BHs per mass bin ordered according to their av- 
erage masses (lines). The asterisks indicate the points where the 
BHs have reached half of their final masses, while empty squares 
represent the time when 90% of the final mass has been accreted. 



tional evidences support the idea of an 'antihierarchical' lu- 
minous activity (star formation or mass accretion onto black 
holes) since it appears to occur earlier in more massive ob- 
jects, a phenomenon generally termed 'downsizing'. 

The evolution of galaxies and BHs are closely con- 
nected, as indicated by the BH-bulge relations fSection l3.2p . 
Therefore, the presence of downsizing in the growth of BHs 
can be understood by taking into account the fact that the 
star formation activity in massive galaxies peaks at higher 
redshifts {z ^ 1.7) than in less massive objects {z ~ 0.5, cf. 
Fig. HI). This phenomenon is reinforced by the effect of AGN 
feedback in masive systems that stops the star formation ac- 
tivity earlier, and consequently the fueling of cold gas that 
feeds the BHs. Downsizing has also been reported in previ- 
ous semi-analy tic mode ls in the study of bo th the growth of 
BH (iBower et al. 2006 .lMalbon et al.|[2007l ') and of elliptical 
galaxies ( Cattaneo et al. 20081 ) . 



3.2 'BH-bulge' relations 

We require our model to reproduce as many observational 
constraints on central BHs as possible. We start by studying 
the relations between the black hole mass and the host bulge 
properties: the 'BH-bulge' relations. 

The first comparison, and the most direct one, corre- 
sponds to the correlation between BH mass and host bulge 
mass. In Fig. U] we co mpare our model re sults (grey points) 
with observations by H aring fc Rixl (|2004 ) , which are shown 
as triangles, along with their best-fitting power law. In order 
to obtain a good agreement with this observational relation, 
we vary the disc instability parameter, edisc, and the param- 
eter /bh (Eq.[5]and Eq. |6] respectively). Even though black 
holes and bulges grow via the same processes, a well ad- 
justed correlation cannot be ensured without a fine tuning 
of these parameters since the BH growth not only depends 
on the amount of available cold gas, as the bulge growth 
does, but also on the BH mass and the virial velocity of the 
halo. 

The relations between BH mass and bulge velocity dis- 
persion, (TBuigc (Fig. [SJ, and between BH mass and bulge 
luminosity (Fig. |6]) also show excellent agreement with ob- 
servations. They are naturally obtained from the choice of 
model parameters motivated by the BH-Bulge mass rela- 
tion. This agreement is not a lways reached by oth er semi- 
analytic models; for instance, iMalbon et al] \2m% show a 
shallower relation between BH mass and bulge velocity dis- 
persion than the observed trend at high values of velocity 
dispersion, while our model produces a good match over a 
larger dynamical range. 



3.3 QSO Luminosity Function 

Whenever a black hole accretes material there is an associ- 
ated bolometric luminosity and, therefore, an active galactic 
nucleus. Several clues on the nature of AGN are provided by 
the study of the evolution of the QSO luminosity function. 

In order to compare the QSO luminosity function ob- 
tained from our model with observational results, we esti- 
mate a QSO absolute magnitude in the optical band Mbj 
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Figure 4. Relation between the BH mass and the host bulge 
mass. Model result s are shown in grey dots. They are compared 
to observations by iHarine &: Kbi j2004l), which are represented 
by triangles with their corresponding error bars; the solid line 
corresponds to the best-fitting power law to these data. 



Bulge 




-2Z -20 -18 -16 -14 

Mb Bulge 



Figure 6. Relation between the BH mass and bulge B— and 
K—hand absolute magnitudes. Model results are represe nted by 
grey points. Compared with observations of .Marconi fc HuntI 
1I2OO3I) (black triangles). 
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Figure 5. Relation between the BH mass and host bulge ve- 
locity dispersion, UBulge- Model results are repres ented by grey 
points. They are co mpared with ob servations of iHaring &: Rixl 
POO^ I (triangles) and lTremain3 lioba l (solid line). 



(blue magnitude in Johnson's system), following Eg. 1121 of 
ICroom et all Hooi), 



-2.66 log,o(Lg„s°^) +79.42 



(12) 



where ieoiom bolometric luminosity in units of Watts 

given by the black hole luminosity defined by Eq. |8l We 
then build the QSO LF at differe nt redshifts and com pare 
its evolution with measurements bv lCroom et~aLl l|2004 ) and 



IWoff et al.1 (|2003l ). The former work presents results for the 
QSO LF from the 2dF QSO redshift survey (2QZ) over the 
redshift range 0.4 < z < 2.1, while the latter uses the 
COMBO-17 survey that contains a large sample of faint 
QSOs in th e redshift range 1 .2 < z < 4.8. Note that the 
LF found bv lWoff et all ([200^ is consistent with results ob- 
tained from surveys such as SDSS and 2QZ. 

The QSO luminosity function is shown in Fig. [T] 
The results of our model (solid lin es) are in good a gree- 
ment with th e observational data of IWolf et al.l (|2003l ) and 
ICroom et al.l |2004) . represented by symbols. These results 
are obtained by the inclusion of the Eddington limit (Subsec- 
tion l2.5|l and a fraction of obscured QSOs. The application of 
the Eddington limit helps to recover the observed QSO LF, 
although a constant normalization offset remains at all red- 
shifts. In order to solve this discrepancy, we find important 
to consider that a fraction /obsc of model QSOs are obscured 
by the gas torus assu med to b e present i n the 'unified model' 
of AGN galaxies (e.g. iFanaroff & Rilc v]|l974l: [Bar thcl 198i; 
iMadau. Ghisellini fc FabianI Il994l : IGunn fc Shanks! [199% : 
we adopt /obsc = 0.8. This fraction is consistent with ob - 
servational results l|Gunn fc Shanksll 19991 : iLacv et al]|2007l ) 
which show that a large number (possibly the majority) of 
QSOs are obscured. 

Model results with /obsc = 0.8 and no Eddington Limit 
are represented as grey dotted lines, showing that this limit 
pro duces changes up to a factor of ~ 6 in the QSO luminos- 
ity. iMalbon et al] (|2007l ) also find that the Eddington limit 
contributes to obtain a good agreement with observations. 
We remark that our simulation box size is only mid-size and, 
therefore, does not allow to reproduce the entire dynamic 
range of observed QSO number densities. 

We stress the fact that we are using a fixed value 
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Figure 7. M odel QSO LF compare d with obse r vation al data by 
ICroom et all 12004) for z <2, and IWolf et~all |20o3) at ^ 3. 
The grey solid lines show the final model including Eddington 
limit with the parameters indicated in Section [2] To obtain good 
agreement with observations it was necessary to assume a frac- 
tion, /obsc=0.8 of obscured QSOs. The resulting model QSO LF 
with no Eddington Limit is represented by grey dotted lines. 
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Figure 8. Upper panel: 6j band galaxy LF. Lower-panel: K- 
band galaxy LF. Solid lines represent results from our new model, 
SAG2, while dashed lines correspond to the previous version, 
SAGl; different sym b ols identify observa tion al LF estimates b y 
iNorbere et"all ||2002D. [Huang et al.l ||2003D andlCole etHI 1I2OOII). 



of /obsc at all redshifts; therefore, it can be said that 
QSO luminosities in our model depend exclusively on the 
physical properties of the AGN phenomenon. We notice, 
however, that at low redshifts the model shows slightly 
higher number densities than the data, possibly imply- 
ing a larger obscured fraction as the redshift decreases. 
There is some observational evidence in this direction, where 
the fraction of obscured QSOs becomes higher at z < 1 
(e.g. iFranceschini. Braito fc Fadda1l2002l : iGandhi fc FabianI 
|2003|). However, this redshift depen dence is not always found 
(e.g. lUeda et al.1 l2003l : iGillil |2004) . being necessary larger 
samples to underst and this tre n d. On the other hand, the 
results obtained by lUeda et al.l (|2003t ) imply that the frac- 
tion of absorbed AGN decreases at high luminosities. Our 
model shows the opposite behaviour; a larger fraction of 
obscured AGN is required at the bright-end of the LF to 
further improve the agreement with the observational data. 

As we have mentioned in Section 12.4.21 the shape of 
the QSO LF is sensitive to the parameter kagn, responsi- 
ble for the regulation of the gas accretion rate onto BHs via 
gas cooling processes. Adopting larger values of kagn leads 
to a higher QSO number density. This increment occurs at 
high redshifts [z > 1) for low and intermediate luminosi- 
ties, while at lower redshifts it mainly affects the bright-end 
of the QSO LF. The latter arises as a consequence of the 
increased impact of gas cooling processes on the growth of 
massive BHs at low redshifts. It is important to remark that 
we do not make assumptions about the QSO lifetime since 
we consider the BH luminosity as the instantaneous rate of 
energy ejected from the BH. 



3.4 Galaxy properties: effects of AGN feedback 

Semi-analytic models are calibrated to reproduce, as best as 
possible, properties of galaxies in the local Universe such as 
the luminosity function, the TF relation, colour-magnitude 
diagram and gas fractions. The luminosity functions in dif- 
ferent colour bands obtained from SAGl are characterized 
by an excess of bright galaxies with respect to the ob- 
served luminosity functions. The TF relation is found to 
lie within the observational er ror determinatio ns, following 
the best-fitting estimation by iGiovanelli et al. I 11997). The 
colour-magnitude relation for early type galaxies managed 
a very good agreement wi th the observational fitting by 
iBower. Lucey fc Ellij l|l992l ). Gas fractions were also found 
to be in very good agreeme nt with observational results by 
iMcGaugh fc de Blokl l| 19971 ). 

In this section, we analyze the advantages of the new 
model SAG2 which now includes AGN feedback, evaluating 
the influence of this process on different galaxy properties, 
and comparing with the results provided by the original ver- 
sion SAGl. 



3.4.1 Luminosity and stellar mass functions 

Luminosity functions in the 6j- and AT-bands are shown for 
SAG2 galaxies in Fig. [S] It can be seen that in both cases 
the model closely follows the trend den oted by the observa- 
tional results from Norberg et al. (|2002|) fo r the &j-band, and 
iHuang et af] l|2003l ) and lCole et al.l7200ll ) for the iC-band. 
At the bright end of the LF, the agreement with observa- 
tions has improved significantly over the previous version 
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Figure 9. SMF for the model galaxies in SAG2 at redshifts 
2 = 4.5, 2, 1 and 0.5 (black solid lines); for comparison, all pan- 
els show the z=0 model SM F (grey soli d lines). Symb ols repre- 
sent o bservational data from lProrv et al.i (,2004 ) and .Drorv et al.l 
1I2OO5I') . 
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Figure 10. Fractions of SAG2 galaxies of different morpholog- 
ical types at 2 = as a function of stellar masses. Model re- 
sults (lines) are compared with observations taken from'Conse licel 
((2006) (symbols): irregular galaxies are represented by dashed- 
dotted lines and triangles, spirals by solid lines and filled circles, 
and ellipticals by dashed lines and open squares. 



of the model, SAGl. This improvement comes from a more 
vigorous quenching of cooling of hot gas at low redshifts 
operating in SAG2, which prevents the presence of very lu- 
minous massive galaxies. 

The masses of model galaxies are able to reproduce the 
observed evolution of the stellar mass function (SMF). The 
dependence of the SMF with redshift is shown in Fig. |9l 
Here we compa re the SMF f r om o ur m odel SAG2 w i th ob - 
servations from iDrorv et al.1 (|2004l ) and iDrory et al.l (|2005| ) 
in four redshift bins, from z ~ 0.5 to 4.5; model results are 
represented by black solid lines. The model SMF is in good 
agreement with observations at all redshifts. The grey lines 
show the z — model SMF, and help to illustrate the lack 
of significant evolution in the SMF from z ~ 0.5. 

Note that our model uses a S alpeter IMF; certain mod- 
els consider different IMFs (e.g. iBaugh et al.l |2004| use a 
top-heavy IMF) in order to reprodu ce the observationa l re- 
sults. It is important to remark that iDrorv et all (|2004 ) and 
iDrorv et al.l ([2005') also assumed a Salpeter IMF to apply 
photometric techniques to determine redshifts and stellar 
masses of the best-fitting templates of galaxy spectra. They 
also analyse two independent fields, the FORS deep field 
and the GOODS-S field, implying that the data take into 
account cosmic variance to some degree and are therefore 
representative of the observed universe at low and high red- 
shifts. 

3.4-S Galaxy morphology 

We analyse the morphology of SAG2 galaxies at z = 0. In or- 
der to determine a fraction of morphological types from the 
model, we use the ratio between the bulge mass and the total 
stellar mass, r=MBuigc/M^, as proposed by Bertone, De Lu- 
cia & Thomas (2007) for this type of analysis. We adopt 
a threshold between ellipticals and spirals rthrcs = 0.95, 



classifying as ellipticals all bulge-dominated systems with 
r > rthrcs, and as spirals those with < r < rthrcs; sys- 
tems with no bulge (r=0) are classified as irregulars. The 
threshold between el lipticals and spirals i s not firmly estab- 
lished; for instance, iBertone et al.l |2003) use rthrcs = 0.7. 
Nonetheless, a higher threshold, as the one adopted here, 
defines very well the elliptical population, allowing them to 
have at the most a stellar disc of 5 per cent the total mass. 

Fig. [TU] shows the fractions of morphological types as 
a function of stellar mass for our model (represent ed by 
lines) and the observational data by IConselicd l|2006l ). rep- 
resented by symbols with errorbars. The data include a 
sample of more than 22, 000 gala:xies ai z < 0.05 with a 
visual morphological cla ssification from the RC3 catalogue 
(|de Vaucourleurg|l991bl ). We can appreciate that the frac- 
tions of different morphological types in the model agree 
very well with the observations over the whole stellar mass 
range. SAG2 shows a clear improvement over the previous 
version of our model, SAGl, in which ~ 80 per cent of mas- 
sive galaxies evolved to be spirals. Starbursts triggered by 
disc instabilities are a key ingredient to achieve a good agree- 
ment with the observed morphological types; disc structures 
are destroyed during these events thus favouring the forma- 
tion of early type galaxies. When including disc instabilities 
in SAGl the percentage of high stellar mass spirals is re- 
duced to ~ 40 per cent. 

Changes in rthrcs, which we varied from 0.7 and 0.95, 
produce a shift in the height of the spiral fraction peak from 
~ 0.55 to ~ 0.75. Spiral galaxies at intermediate stellar 
masses are the more sensitive population to rthrcs. Ellip- 
tical galaxies are the dominant type at high stellar masses 
regardless of the value of this parameter. 

In order to test the reliability of the morphological clas- 
sification described above, we apply the criteria given by 
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iDe Lucia et~all l|2004l ). in which spiral galaxies satisfy the 
condition 1.2 < Mb, Bulge — MB,Gai < 2.4; we also con- 
sider t he definition of Hubble types given by 'Sprin gel et al.l 
Our conclusions are not significantly affected in any 
of these cases. 

Recent works on the analysis of morphological types 
in semi-analytic ga l axies obtain slightly different results. 

Lucia fc Blaizod l|200'it fail to predict t he presence of 
high m ass spiral galaxies. On the other hand. iBertone et al.l 
(|2007f ) show a good agreement with observations at all stel- 
lar masses although their model produces a significant excess 
of bright galaxies in the LF and a smooth colour distribution 
with most of their galaxies residing in a red sequence. The 
morphological aspect of model galaxies is quite sensitive to 
the different parametrizations of the processes affecting the 
star formation, including the effects of disc instabilities and 
major and minor mergers. Thus, it is not surprising that 
different semi-analytic codes lead to different dependencies 
of morphological type fractions with stellar mass. 



3.4-3 Galaxy colours 

The AGN feedback has an important effect on galaxy 
colours. Fig. [11] presents B — V colours as a function of 
stellar mass for the models SAGl (left panel) and SAG2 
(middle panel), where galaxies are separated in ellipticals 
(red symbols) and spirals (blue symbols) following the crite- 
ria outlined in Subsection 13.4.21 The right panel shows the 
same SAG2 galaxy population splitted into central (black 
symbols) and satellite galaxies (green symbols). 

SAG2 galaxies show that the inclusion of AGN feed- 
back produces galaxies with B — V colours qualitatively 
compa tible with measured colour distributions (see for in- 
stance iBaldrv et al.l I2OO6I ) including a bimodal behaviour 
and a negligible population of massive (> lO"/i-^M0)blue 
galaxies. This is a clear difference with respect to SAGl 
where there is a population of very massive, blue spiral 
galaxies. As it was mentioned in Subsection 13.4.21 this is 
a consequence of the generation of new spheroids by the 
action of the disc instability process included in SAG2. 

The right panel of Fig. [11] indicates that the red branch 
is mainly composed by satellite galaxies, characterized by 
the fact that their hot gas reservoir has been stripped when 
they were accreted by the central halo, aging quiescently 
until merging with a central galaxy. Massive central galaxies 
also become red as a result of a combination of different 
effects. On the one hand, the reduction or suppression of 
gas cooling leaves a galaxy with a limited cold gas reservoir 
that is eventually exhausted; it then undergoes a passive 
evolution unless a merger with another galaxy takes place. 
However, the possibility that this occurs decreases at low 
redshifts, as it is shown in Fig. [1] On the other hand, if 
a merger occurs, it is very unlikely that satellite galaxies 
will contribute enough cold gas to produce a starburst that 
will appreciably change the colour and mass content of the 
central galaxy, since the gas cooling process does not occur 
in satellites. Hence, mergers would more likely only increase 
the mass of the stellar population without giving place to 
an important generation of new stars. 



4 FURTHER IMPROVEMENTS TO THE 
SEMI-ANALYTIC MODEL OF GALAXY 
FORMATION 

In this section, we explore the effects of two new ingredients 
that can be added to the model. We first take into account 
an additional form of feedback from AGN, associated with 
the black hole growth during starbursts. Then, we focus on 
the analysis of the effects of different prescriptions for the 
dynamical friction time-scales affecting galaxy mergers. 



4.1 AGN Feedback in starburst events 

The importance of AGN feedback in galaxy evolution is un- 
deniable, as theoretical and observational works indicate, 
although this is still a poorly understood phenomenon and 
is at present a matter of debate. It is mostly accepted that 
the quiescent gas accretion produces feedback from the ac- 
tive nucleus, but there is no general consensus on the fate 
of the energy produced by the BH mass accretion during 
starbursts. Therefore, in this section, we use our new semi- 
analytic model SAG2 to explore the effects of AGN feedback 
during violent events (i.e., starbursts) on galaxy properties. 

There are observationa l clues about the infl uence of 
AGN feedback in starbursts. [Donahue et al.l l|2005h propose 
a dichotomy in galaxy clusters between those which are radio 
active and show strong central temperature gradients (indi- 
cating the action of AGN feedback), and radio quiet with no 
temperature gradients. Nonetheless. iGastadello et al.l (|2007l ) 
find a counter-example in the cluster AWM 4 characterized 
by the absence of a temperature gradient, as indicated by 
the XMM data, and the presence of an active radio source 
at 1.4 GHz. The combination of these features cannot be ex- 
plained if AGN activity is regulated only by feedback from 
gas cooling processes. 

We explore the effects of AGN feedback during star- 
bursts following two different possibilities: (i) the outflowing 
material is directly transferred to the hot gas surrounding 
galaxy, and (ii) the energy produced is high enough to ex- 
pel material out from the galaxy halo without affecting any 
baryonic component associated to the galaxy. We refer to 
these two mechanisms as retention and ejection modes, re- 
spectively; the former has been used in both SNe feedback 
and AGN feedback in radio mode throughout this work. 
The absence of definitive evidence favoring one scenario or 
the other only allows us to analyse and compare the out- 
comes from each option. Note that the ejection mode does 
not operate exactly i n the same way as in, for example, 
iDe Lucia et al.l (|2004[ ). where the ejected gas is eventually 
reincorporated to the dark halo. AGN feedback in starburst 
mode heats the cold gas already present in the galaxy with 
a rate 2 x Z/Bn/Vvir- The black hole luminosity, -Lbh, given 
by Eq. [9] now involves the mass accreted in the starburst 
mode, as given by Eq. 1101 

Fig. ll2l shows the relation between black hole and bulge 
masses considering AGN feedback produced when the ac- 
cretion occurs as a result of starburst events, for the case 
of the ejection mode. The results from the retention mode 
are not shown in this figure since the general trend does 
not vary appreciably with respect to the ejection scenario, 
and only produces a larger scatter of BH masses. As can be 
seen, AGN feedback during starbursts produces only slight 
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Figure 11. Colour distribution for ellipticals (red symbols) and spirals (blue symbols) type galaxies obtained from models SAGl (left 
panel), and SAG2 (middle panel). SAG2 galaxies separated into centrals and satellites are also shown (right panel, black and green 
symbols, respectively). 




Figure 12. Relation between BH and bulge masses for model re- 
sults including AGN feedback during starbursts in the ejection 
mode (grey points) compared with observational estimates by 
iHaring fc Rb3 ||2004|) (triangles). 



changes with respect to the model where only radio mode 
AGN feedback is applied (see Fig.|4]). The main difference is 
a shallower relation produced by the complementary effect 
of AGN feedback in starbursts, that reduces even more the 
available cold gas for bulge growth so that, for small BIfs, 
bulge masses are smaller. Note that this does not affect BH 
masses as significantly since BHs grow mainly through merg- 
ers and disc instabilities (cf. Fig. In any case, the poor 
statistics in the observational relation does not allow us to 
determine which of the three versions of our model (AGN 
feedback from gas cooling, or AGN feedback during star- 
bursts for the retention/ejection mode) provides the best 
agreement. 

We have to take into account that AGN feedback in 
starburst mode is an additional ingredient that changes the 



mass and metallicity of the hot gas, which in turn affects 
the cooling rate and thus the stellar formation rates. Fig. 1131 
shows the global cosmic star formation rate as a function 
of redshift for the models SAGl (dashed Une), SAG2 (solid 
line) and SAG2 including AGN feedback in starbursts for 
the retention (dotted line) and ejection (dot-dashed line) 
modes. The latter two models use the same set of parame- 
ters as SAG2. Comparing SAG2 with SAGl, we can appre- 
ciate that, for z > 3, the SFR is larger in SAG2 as a result 
of the disc instability and starbursts in minor mergers in- 
cluded in this model. The feedback in radio mode becomes 
progressively more effective for lower redshifts, increasing 
the reheated mass (solid grey line in lower panel of Fig. I13|l 
and appreciably reducing the SFR. When including AGN 
feedback in starbursts in the ejection mode, the SFR is re- 
duced at all redshift both with respect to SAGl and SAG2, 
keeping the similar redshift dependence as in SAG2. This is 
a consequence of the removal of gas by the ejection mode, 
which is no longer available to cool and form stars. In the 
case of the retention mode for the AGN feedback in star- 
bursts, the gas is mainly reheated by the radio mode, as 
shown in the lower panel of Fig. 1131 This is due to a larger 
reservoir of hot gas available for gas cooling due to the re- 
tention of the gas heated by the AGN. This produces a much 
higher SFR in this mode with respect to the other models 
shown in the figure. 

Three models, SAGl, SAG2, and SAG2 with retention 
mode AGN feedback in starbursts, show a re asonable agree- 
ment with the observational data compiled bv lHopkins et all 
(j2006i ) (grey symbols with error bars) at all redshifts. How- 
ever, SAGl results show a too high SFR at z = 0. On 
the other hand, SAG2 with AGN feedback in starbursts 
in the ejection mode, shows a much lower star formation 
history (SFH) at all redshifts, although it should be noticed 
that some observational data points are consistent with even 
lower values of SFR at z > 2. Note that SAG2 provides a 
good agreement with the oberved SFH at all redshifts. Only 
at ~ 1 the model shows slightly lower values than the dis- 
tribution of observational results. How ever, recent studies by 
IWilkins. Trentham fc Hopkins! (|2008l ) show that direct mea- 
surements of star formation rates (e.g. iHopkins et al.|[2006l ) 
seem to be overestimated with respect to the SFR values ob- 
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Figure 13. Upper panel: Cosmic star formation rate for SAG2 
(solid line), SAG2 including AGN feedback during starburst 
(with the same parametrization as SAG2; dotted and dot- 
dashed lines for the retention and ejection modes, respectively) 
an d SAGl (dash e d line ), compared to observations compiled 
by iHopkins et al] ||2006| ) (grey triangles with errorbars). Lower- 
panel: Reheated mass from the different feedback mechanisms as 
a function of redshift for the different AGN feedback models. The 
grey lines represent the radio mode and the black lines the star- 
burst mode. 
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Figure 14. bj- (upper panel) and X-band (lower panel) lumi- 
nosity functions for three models: when AGN feedback is acting 
only during gas cooling process (solid line), when AGN feedback 
is also produced in starburst events with ejection of the reheated 
material (dotted line); and when the material is retained in the 
dark halo in which the galaxy resides (dashed line). In the ejec- 
tion mode we use the same set of parameters as in SAG2, while 
a value of kagn = 2 x 10"'^ Mq yr~^ is adopted in the retention 
mode. 



tained from stellar mass functions. Taking this into account, 
the agreement between our model SAG2 and observations 
at 2 ~ 1 would be even more consistent. 

We also analyse how the inclusion of AGN feedback 
in starburst mode affects the galaxy luminosity function. 
Fig. [14] shows the 6j- (top panel) and ii'-band (bottom 
panel) LPs for the three different models: SAG2 (solid line), 
and SAG2 including feedback in starbursts in the retention 
(dashed line) and ejection modes (dotted line). As can be 
seen in Fig. [TSl S AG2 gives lower SFRs when AGN feedback 
in starbursts with ejection is included. Therefore, in this case 
the LF is slightly shifted towards fainter luminosities by up 
to half a magnitude. As the agreement with observed results 
is not badly affected by this, we do not need to recalibrate 
the model parameters in this case. In the retention mode, the 
radio mode feedback efficiency, the kagn parameter (Eq.O, 
needs to be increased in order to reproduce the observed 
bright-end of the LFs; we adopt kagn — 2 x 10~^ M© yr~^. 
After applying this recalibration, the retention mode pro- 
duces a good agreement with observations. 

The morphology of model galaxies (estimated as in Sec- 
tion I3.4.2|l is affected by the AGN feedback during star- 
bursts. The ejection mode improves the agreement with 
the observed m orphological fractions at low stellar masses 
(|Conselicell2006l ). producing a non-zero spiral population at 
this mass range. This mode also mantains the good agree- 
ment at intermediate and high stellar masses shown by 



SAG2. On the other hand, the retention mode fails to re- 
produce the observed morphological fractions at all masses, 
even after the feedback efhcency, kagn , is increased in order 
to match the observed galaxy LF. 

The failure of the retention mode can be used as an in- 
dication that, if present, AGN feedback during starbursts 
should be strong enough to expell the reheated material 
away from the galaxy halo. The only drawback to this ejec- 
tion model is that it may produce slightly low values for 
the SFH with respect to the observational trend, which in 
any case still shows data points that lie below the model 
prediction at high, z > 2, redshifts. 

The implementation of AGN feedback in starbursts 
is still crude since it considers the same efficiency as in 
the AGN feedback acting during gas cooling. It is possible 
though, that the physical considerations needed for these 
two different accretion processes, via gas cooling and during 
starbursts, are not the same. For instance, instead of using 
different BHs feeding mechanisms to distinguish bet ween dif- 
ferent feedback modes, the model described bv iSiiacki et al.l 
(|2003) uses a threshold value for BH accretion rate above 
of which the feedback is similar to our starburst mode al- 
though with a lower efficiency. We leave the detailed study 
of AGN feedback during starbursts to a forthcoming paper. 
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Figure 15. fej- and if-bands galaxy luminosity function for three 
di fferent merger ti me-sc ales; the results using the prescription 
by iLace v fc Co"i3 l ll993t) are shown by solid lines, those from 
Ijiang et al. I | |2007^ and iBovlan-Kolchin et al.l fcOOTt) are repre- 
sented by dashed and dotted lines, respectively. 



merger frequencies. Both major and minor mergers are af- 
fected, although the changes are more important for the mi- 
nor mergers. This translates into a shift in the period of max- 
imum merger activity towards redshifts closer to z ~ 1.5, 
much lower than that obtained for LC93, which occurs at 
z ~ 2.5. In turn, this affects the disc instability frequency 
presented in Fig. [T] producing an order of magnitude in- 
crease in the number of events. This is a result of the longer 
time-scales over which galaxies can form stars and accrete 
gas into the disc. The higher number of disc instabilities 
and the delayed action of mergers compensate each other, 
and no recalibration of the model parameters is needed in 
order to obtain a good agreement between model and ob- 
servational galaxy properties such as the z = galaxy LF 
and bimodal galaxy colour distributions, stellar masses at 
different redshifts, 'BH-bulge' relations, and the SFH. 

The new time-scales affect galaxy morphologies produc- 
ing a 10 per cent excess of elliptical galaxies at int ermediate 
stella r masses with respect to observations (e.g. IConselic3 
l2006h . A possible reason for this slight excess of spheroids 
is that the frequency of disc instability events increases pro- 
ducing a somewhat larger population of smaller elliptical 
galaxies than SAG2 and the observations. At the high mass 
end, the J07 prescription produces a slightly better agree- 
ment with the observed morphological fractions. 

As it was mentioned above, the effect of the new merger 
time-scales is difficult to detect in the galaxy statistics stud- 
ied in this paper. However, the study of quantities that are 
affected by the detailed stellar mass growth may show more 
important differences, as would probably be the case for the 
amount of stellar mass formed in accreted satellites. 



4.2 New prescriptions for merger time-scales 

All the results presented in this work have been obtained 
considering that galaxy mergers are driven by the dynam- 
ical friction process, for which we have adopted the pre- 
scription given by Lacey & Cole (1993, hereafter LC93). We 
now take into account the new dy nami cal friction time-scales 
propo sed by IJiang et al.l (|2007l ) and IBovI an-Kolchin et al.l 
l|2007l ) (hereafter J07 and B07, respectively). Both stud- 
ies use numerical simulations to measure merger time- 
scales and compare them to theoretical predictions based 
on the Chandrasekhar form u la. Th ey find that, in contrast 
to iNavarro. Frenk fc White! (|l997l ) , the theoretical merger 
time-scales are not in agreement with the simulated values. 
Equations (5) from J07 and (5) from B07 provide new fit- 
tings to the merger time-scales that we test in our model. 
The difference between these two fittings resides in the 
weaker dependence on circularity in the prescription of J07 
compared with the one by B07. The overall conclusions from 
both studies are the same, that the LC93 merger time-scales 
are underestimated for minor mergers and overestimated for 
major mergers. 

The effects of the three different merger time-scales on 
the galaxy population can be readily seen in Fig. 1151 where 
we show the luminosity function in the bj- and K-hands (up- 
per and lower panels, respectively; using the SAG2 param- 
eter set) for the three different approximations. The three 
models are in good agreement with observations in both 
bands. 

Changing the merger time-scales according to the new 
prescriptions of J07 and B07 produces a decrease in the 



5 SUMMARY AND CONCLUSIONS 

As several recent observational results indicate, AGN feed- 
back has become a very important element in our un- 
derstanding on galaxy evolution (e.g. iNulsen et al.l |2007| . 
[Schawinski ot aL 120071 ). In this paper, we study the effects 



of AGN feedback on the formation and evolution of galaxies 
by using a combination of a cosmological A'^body simulation 
of the concordance ACDM cosmology and a semi-analytic 
model of galaxy formation. We have developed an AGN 
feedback model based on different im plementations con- 
sidered in recent semi -analytic models (jCroton et al.l |2006| . 
iBower et al]|2006l and lMalbon et al.ll2007l '). which are based 
on both, observations and a phenomenological picture of the 
physics governing the extremely massive, central BHs, and 
of the gas ejected by these central engines back to the inter- 
stellar or even the intergalactic medium. Our semi-an alytic 
model is an improvement of the one described by ICoral 
(|2006D which involves gas cooling, star formation, galaxy 
mergers and SN feedback, with the the advantage of in- 
cluding a proper treatment of the chemical enrichment of 
baryons. The new version of this model (SAG2) benefits 
from the inclusion of AGN feedback. In our model, the AGN 
activity is triggered by the accretion of cold gas onto BHs, 
which grow via three distinct channels: (i) mergers between 
two pre-existing black holes, (ii) mass accretion during gas 
cooling processes, and (iii) gas accretion during starbursts. 
The latter are produced during major and minor mergers, 
and disc instabilities, being the last two processes new in- 
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gredients to our model. Black holes produce a luminosity 
which is proportional to the gas accretion rate. When the 
accretion is produced during the gas cooling process, the 
resulting AGN produces feedback (radio mode feedback). 
It is also considered the possibility that AGN generated by 
the BH growth during starbursts also produce feedback. An- 
other novel feature analysed in our model is the influence of 
new prescriptions for merger time-scales on galaxy proper- 
ties. 

We summarise the main results of this work. 



• The growth history of galaxies is dominated by both 
major mergers and disc instabilities, while minor mergers 
play only a minor role (Fig. The frequency of starbursts 
peaks between z ^ 2 — Z when produced by mergers, and 
between 2 ~ 1 when produced by disc instabilities. 

• The accretion history of black holes shows small varia- 
tions between high and low mass BHs. For high mass BHs 
at z = 0, M-BYi{z = 0) > 10** Mq, the accretion is dominated 
by disc instabilities at all redshifts. Instead, mergers dom- 
inate at 2; > 2 for lower masses, with a larger influence of 
disc instabilites at lower redshifts (Fig.[21). Regardless of the 
BH mass, the accretion via gas cooling processes is negligi- 
ble. Higher mass black holes assemble earlier in the history 
of the Univer se (Fig. El, in agreeme nt with the downs izing 
scenario (e.g. iMarconiet ah 2004, S hankar et al.ll2004 l. 

• We simultaneously obtain good agreement with ob- 
servations in several 'BH-bulge' relations: BH-bulge mass 
(Fig. 131 , BH mass as a function of bulge velocity dispersion 
(Fig. [5]) and bulge luminosity (Fig.[6ll; these agreements were 
achieved using a fixed value of the parameter /bh (Eq. [6]) 
tuned to match the observed BH-bulge mass relation. 

• Our model predicts that a large fraction (80%) of 
model QSOs need to be obscured in order to reproduce 
the observed QSO LF (Fig. [Tjl. Many studies have indi- 
cated this fac t but fail to provide ac c urate estimates o f this 
fraction (e.g. iGunn fc Shanks! Il999l . iLacv et al.1 l200'iD . At 
low redshifts the model shows slightly higher QSO num- 
ber densities than the data, indicating the possibility that 
the obscured fraction needs to become larger as the red- 
shift decreases, as suggested by some observational results 
fe.g. iFranceschini. Braito fc Fadda1l2002l : [Gandhi fc FabianI 
I2OO3I V 

• SAG2 reproduces many of the observed properties of 
the local galaxy population that were also reproduced by 
the previous version of the model (e.g. the TuUy-Fisher re- 
lation, and the relation between gas fraction and absolute 
Mb magnitudes) . The observed bright-end of the galaxy lu- 
minosity functions in the 6j- and A"-bands (Fig. [5} have 
been successfully recovered by SAG2 galaxies, being one of 
the main successes of the model. Our model also agrees very 
well with observational determinations of the stellar mass 
function out to very high redshifts [z ~ 4.5, Fig. [9]), as well 
as the global star formation rate history (Fig. llSp . The latter 
clearly demonstrates how AGN feedback efficiently reduces 
star formation by quenching the gas cooling process at low 
redshifts. 

• Our model re produces the ab undances of observed mor- 
phological types (|Conselicel |2006| ) at all stellar masses. We 
notice that the morphological fractions are very sensitive to 
the AGN efficiency, which was set by comparing model and 



observed galaxy LFs, an independent statistic of the galaxy 
population. 

• Our model galaxies show a bimodal B- V colour distri- 
bution up to stellar masses « lO^^'^'^^M©. At higher 
masses, the distribution becomes unimodal, presenting al- 
most only red, elliptical galaxies. This constitutes one of 
the main effects of AGN feedback since the original model 
shows a bimodal distribution for all stellar masses with an 
important massive, blue spiral population. The transition 
to a unimodal distribution is due to the passive evolution of 
galaxies produced after exhausting their cold gas. 

• We have implemented two different channels of star- 
burst feedback: 'retention mode' and 'ejection mode' (in- 
dicating that gas is either retained or totally ejected from 
the galaxy halo, respectively). We found that the ejection 
mode provides a good agreement with all the observational 
constraints mentioned above. This mode does not require a 
recalibration of model parameters with respect to those used 
in SAG2, in contrast to the retention mode which demands 
a higher AGN feedback efficiency (kagn) in order to repro- 
duce the observed galaxy LF; however, even after recalibrat- 
ing the model, the morphological fractions in the retention 
mode disagree with the observational measurements. This 
indicates that AGN feedback during starbursts needs to be 
strong enough to expell the reheated gas from the galaxy 
halo. 

• We implemented the ne w prescripti ons for merger 
time-scales presente d by Ijiang et al.l l|2007h and 
iBovlan-Kolchin et al.l |2003) These prescriptions pre- 
dict that major merger time-scales are signifi cantly shorter 
than those arising from the estimation of iLacev fc Cold 
(|l993l ). As result of this, galaxies have more time to accrete 
gas and stars in the disc, thus producing a significantly 
higher number of disc intability processes than what is 
shown in Fig. [T] However, the net effect from applying 
the new time-scales is difficult to detect, and the z = Q 
galaxy properties (e.g. galaxy LF, Fig. llSp remain almost 
unaffected. 

Previous semi-analytic models were only focused on 
the analysis of par ticular properties of galaxies and BHs. 
ICroton et al.l (I2OO6I') pre s ent re s ults on galaxy p r opert ies at 



z=0, while 



Bower et al.1 (|2006l ). ICattaneo etalj l|2006l ) and 



iMenci et alT i 20061 ) also take into account properties at high 



redshift but without a compl ete analysis of black hole be- 
haviour. On the other hand, iMalbon et al.l (|2007h make a 
complete analysis of BHs but do not present an analysis 
on galaxy properties. The combination of the different pre- 
scriptions of the black hole growth and the associated AGN 
feedback taken from these theoretical works allowed us to 
develop a model that is able to reproduce several observa- 
tional constraints simultaneously, focusing our analysis on 
different important results and predictions at low and high 
redshifts. We obtain similar results for the galaxy LF and 
colour distributions as previous models, but also obtain im- 
portant differences in the evolution of the black hole accre- 
tion rate; these differences may prove to be useful to obtain 
more clues on the appropriate modelling and parametriza- 
tion of starburst processes, chemical enrichment and AGN 
feedback models. 

We have presented a new model of galaxy formation and 
evolution where we consider different mechanisms for the 
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growth of black holes and galaxies, embedded in a ACDM 
cosmology. Among the several indicators that follow ob- 
servational constraints, it has emerged that the black hole 
and bulge growth may seem 'anti-hierarchical', in agreement 
with the down-sizing scenario, even though the roots of their 
growth in our model is the hierarchical clustering scenario. 
This phenomenon deserves a detailed analysis which we will 
tackle in a forthcoming paper. Further modifications to the 
model are being developed in order to improve the treat- 
ment of the black hole physics, which involves the BH spin 
and different physic al models of accretion (e.g. lBondi|[l952l . 
Naravan. Gsjcia fc McClintock 1997) and disc warping (e.g. 
Shakura &: Sunvaevlil973l . iKing et al.l i2005V The implemen- 
tation of these processes will give us a powerful tool to study 
into more detail the origin of the obscuration of QSOs and of 
the radio-loud QSO population, the relation between black 
hole properties and galaxy morphology, and several other 
interesting observed AGN features. This will provide more 
links between the diverse AGN phenomenon and new un- 
derlying still unknown processes involved in the formation 
and evolution of galaxies. 
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